Electron-phonon coupling and its evidence in the photo-emission spectra of lead 
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We present a detailed study on the influence of strong electron-phonon coupling to the photo- 
emission spectra of lead. Representing the strong-coupling regime of superconductivity, the spectra 
of lead show characteristic features that demonstrate the correspondence of physical properties in 
the normal and the superconducting state, as predicted by the Eliashberg theory. These features 
appear on an energy scale of a few meV and are accessible for photoemission only by using modern 
spectrometers with high resolution in energy and angle. 

PACS numbers: 71.30.+h 74.25.Jb 74.70.Ad 79.60.-i 



A striking evidence of electron-phonon interaction in 
solids is the existence of superconductivity. Bardeen, 
Cooper, and Schrieffer showed Q that even a weak 
electron-phonon coupling is able to condense two elec- 
trons to a so-called Cooper-pair, which is the basic pre- 
requisite of the superconducting ground state and the key 
feature of the BCS model. The BCS theory describes suc- 
cessfully the superconducting properties of many solids, 
e.g. Al, V, or VaSi, where the electron-phonon coupling 
is sufficiently weak. Other conventional superconductors, 
e.g. Pb, Hg, or NbaGe, show significant quantitative and 
qualitative deviations from the predictions of the BCS 
model . These systems are usually classified as strong- 
coupling superconductors. 

The theoretical approach for the explanation of strong- 
coupling superconductors is based on the Eliashberg 
equations, with the coupling function a 2 F as the central 
property [3j. This so-called Eliashberg function can be 
calculated e.g. by first principle methods from the elec- 
tronic wave functions, the phonon density of states, and 
the electron-phonon coupling between two Bloch states. 
However, a 2 F is a much more universal quantity and de- 
termines also the influence of the electron-phonon inter- 
action on the normal state properties, e.g. the electrical 
resistivity, the electronic heat capacity, and — as a spec- 
troscopic feature — the electron-phonon contribution 
Tei—ph to the intrinsic quasi-particle linewidth, which can 
be determined e.g. by photoemission spectroscopy. In 
the Green's function method, the influence of electron- 
phonon interaction is expressed as a contribution to the 
complex self energy of the conduction electrons S e ;_ p ^, 
which can be calculated from the Eliashberg function 
a 2 F ||. The real part describes the induced band renor- 
malization, where the imaginary part gives the quasi- 
particle line width, equivalent to the reciprocal hole life- 
time r. In particular, at very small energies (i.e. close 
to the Fermi level), the real part of E e ;_ p ^ is linear in 
energy, and A = — 8^.eT, e i- p h(w) / 5uj at u = is usu- 
ally called the mass enhancement factor. In addition, A 
defines the slope of the temperature dependence of the 
quasi-particle linewidth Y e i- p h = 25m£ e ;_ p / 1 w 27rAfcsT 
for temperatures above the Debye temperature Orj. 



Photoemission spectroscopy (PES) is a versatile exper- 
imental method to study the electronic structure of solids 
and has been applied to many high- T„ materials |5j,|6J and 
conventional superconductors 0, H, 0, 0, El . Electron- 
phonon coupling in metallic systems has been studied 
in detail by PES on low-dimensional electronic states at 
surfaces, e.g. Shockley-states 111 [H Q H H [H E3 
or quantum-well states 0, l20| . PES investigations 
on electron-phonon effects in three-dimensional solids, 
however, are rare. One particularly interesting three- 
dimensional model system is Pb, which among the 
conventional superconductors has quite un-conventional 
physical properties. 

In the case of Pb the two most relevant energy scales on 
which the electron-phonon features appear in the spec- 
tra, the Debye energy huirj = &bOd an d the gap width 
Ao at T = 0, define prerequisites that can be met by 
PES experiments today {ftwo = 7.6 meV, A = 1.4 meV). 
The instrumental resolution must be at least of the same 
order of magnitude to give usable spectral information 
on these features, i.e. AE < 3 meV. Furthermore, if one 
wants to measure the band renormalization on single- 
crystalline samples, the angular resolution must be suf- 
ficient to resolve the band dispersion of states close to 
the Fermi level. Since a few years, photoelectron spec- 
trometers fulfil these requirements and PES has become 
a unique method yielding access to electron-phonon fea- 
tures in both the density of states and the fc-resolved 
electronic band structure. 

The photoemission data presented here have been mea- 
sured using He l a radiation from a monochromatized he- 
lium discharge lamp (hv = 21.22 eV) with an energy 
resolution of < 3 meV. The angular resolution was typ- 
ically « 0.3° and the sample temperature could be var- 
ied from room temperature down to about 4.5 K (see 
Ref. 0,0] for details of the experimental setup). The 
samples were polycrystalline and single-crystalline pel- 
lets (Surface Preparation Laboratory, Netherlands), ther- 
mally connected to the sample holder by soldered in- 
dium. Clean single-crystalline surfaces were prepared in 
situ by careful Ar sputtering and subsequent annealing 
at < 150°C. 
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FIG. 2: Temperature dependence of the lifetime width V for 
energies uo > u)r> (circles). The dashed and the solid line rep- 
resent the calculated electron-phonon contribution at Ef = 
and E = fta> ma x = 12 meV, shifted vertically to match the 
experimental data. The inset shows a 2 F(u) as used for the 
calculation (from Ref. |22j|1. 



FIG. 1: Real and imaginary part of the electron-phonon self 
energy E, experimental data (circles) and calculated from Ref. 
[5J (full lines). Upper inset: measured band dispersion close to 
Ef in the normal phase (T«8 K); the dashed line gives the 
result of a linear fit to the data points at higher energies (Eb > 
20 meV). The experimental data points have been determined 
from the maximum position of a Lorentzian fit to the angular 
distribution curves (ADC, see lower inset with ADC at Eb = 
10.5 meV). 



On such a surface, one is able to investigate the elec- 
tronic band structure of states close to the Fermi level by 
angular resolved PES (ARUPS). The inset in the upper 
panel of Fig. Ogives a volume band dispersion Ej- close to 
the Fermi level on a Pb(110) surface in the normal state 
(T = 8 K). With He I Q , this band, forming the locally 
tubular Fermi surface of Pb, can be reached slightly off 
the high-symmetry direction TK [23, or, more precisely, 
at k\\, x = 0.82A" 1 along TK and k\\ >y = O.llA -1 parallel 
to FX. The data points (circles) were determined from 
the position of the peak maxima in the angular distri- 
bution curves (ADC) at the energies close to the Fermi 
level Ef (see lower inset). At higher energies, i.e. sig- 
nificantly higher than the Debye energy Hujd — 7.6 meV, 
the dispersion can be approximated quite well by a linear 
function of the momentum fc, representing the dispersion 



ek with no renormalization. However, below E « 12 meV 
the band shows the well known correction due to electron- 
phonon interaction, which is defined by the real part of 
the self energy 5ReS(i?fe) = Ek — fife- The upper panel 
shows the experimental result (T = 8 K) for the deviation 
from the linear dispersion compared to the calculated 
5ieE from Ref. at T = 11 K. The comparison shows a 
good qualitative agreement, in particular the position of 
the maximum at w 7 meV and the slope at Ep (given by 
the coupling constant A = 1.55) are the same for experi- 
ment and theory, within the experimental uncertainties. 
The difference in maximum height can be explained by 
the finite experimental resolution. Up to now, such a 
renormalized quasi-particle dispersion has been observed 
only for low-dimensional systems, e.g. the Shockley-state 
on Mo(110) 0] or high-T c superconductors 

By an analysis of the ADC linewidth as a function of 
the binding energy one can extract the energy depen- 
dence of the imaginary part SmE of the self energy (see 
e.g. Refs. 0, lisjh The lower panel of Fig. ^ shows 
the experimental energy dependence in comparison to 
SmSfBt) from Ref. H Since Pb is a three-dimensional 
system, the contribution of the final state to the total 
linewidth can not be neglected [3| • This additional con- 
tribution is not related to electron-phonon interaction 
and therefore does not show an energy or temperature 
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FIG. 3: Lower panel: real and imaginary part of the gap 
function, interpolated at T = 4.4 K from the calculated gap 
functions at T = and 0.98T C given in Ref. pfij . Upper panel: 
resulting spectrum, raw normalized density of states N/No 
(solid line) and broadened by a convolution with a Gaussian 
(FWHM A_B = 2.8 meV) to simulate the resolution function. 



S — f> exp. data at T=4.5 K 

model function 

I (difference) 



20 15 10 5 
binding energy [meV] 



E F =0 




FIG. 4: Comparison of the modelled function from Fig. [3] 
(solid line) with the experimental data (circles). The black 
bars at the x-axis indicate the difference between experiment 
and theory. The inset shows a blow-up of the peak and dip 
structure. 



dependence in the investigated range. Thus, it is rea- 
sonable to simply shift the theoretical result by 82 meV 
to match the photoemission data. Apart from this off- 
set and experimental limitations, the agreement between 
experimental result and theory is good: below the maxi- 
mum phonon energy hw max ~ 10 meV (see inset of Fig. 
|2J) there is a continous increase of the linewidth up to 
a net change of w 13 meV. As known from the Debye 
model, above Hu> = 10 meV the electron-phonon contri- 
bution to the linewidth remains constant. 

The photoemission linewidth is also characteristically 
dependent on the temperature. Fig. [21 shows the exper- 
imental linewidth from ADCs for u> > u>d (circles); the 
dashed and the solid line give the theoretical results at 
Ep and Eb — tv^max-, respectively. Obviously, the calcu- 
lated temperature dependence describes the experiment 
quite well. For the whole investigated temperature range, 
the theoretical curve lies within the error bars of the ex- 
perimental values and, in particular, for high temper- 
atures the experimental linewidth follows the predicted 
linear behavior given by r e /_ p ^ = 27tA/cbT ', again with 
A = 1.55. This result confirmes that the contribution of 
the final state to the total photoemission linewidth does 
not considerably depend on the temperature. 

Obviously, there is clear evidence for the strong 
electron-phonon interaction in the photoemission spectra 
of Pb in the normal state. In addition, there are char- 
acteristic spectral features in the superconducting state 



of Pb: a peak-and-dip structure, that has been experi- 
mentally observed first by tunneling techniques [2r3 | and 
reproduced recently in photoemission data ||. Theo- 
retically, the density of states of a strongly coupled su- 
perconductor can be fully described by an energy and 
temperature dependent, comple x gap f unction A(E, T) 
given by the Eliashberg theory [271 l28j . The imaginary 
part of A is the consequence of the damping of the quasi- 
particle excitations caused by the electron-phonon inter- 
action. With this gap function, the spectral density of 
states of a strongly coupled superconductor is described 
by N(E,T)/N = 3?e [e/^/E 2 - A 2 (£, T)}, giving as 

a special case the dip-less BCS density of states when A 
chosen real and constant in E. 

A theoretical gap function for Pb at finite tempera- 
tures can be found in Ref. 25, which gives the energy 
dependence for A(E,T) explicitly for two temperatures 
Ti = and T 2 = 0.98T C . To obtain the gap function at 
the temperature of the present experiment (T = 4.4 K), 
we simply interpolated between these two pairs of curves 
A(E,T) = (l-a(r))A(E,T 1 )-r-a(T)A(E ) T 2 ), consid- 
ering the temperature dependence of the BCS gap [2j|. 

We model the photoemission spectrum by convoluting 
the result with a Gaussian to describe the broadening 
from the finite experimental resolution. The interpolated 
gap function and the resulting spectrum for T — 4.4 K 
are displayed in Fig. [31 

Fig. 01 shows a comparison of the model function with 
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the experimental data at T — AA K, normalized to unity 
at E = 27 meV where the spectrum becomes flat; the 
inset shows a blow-up of the range of the dip. The solid 
line represents the model function from Fig. |3 the bars 
at the x-axis represent the difference to the normalized 
experimental spectrum. The agreement between exper- 
iment and model is striking. Without using any free fit 
parameter, all important features — the position of the 
gap edge, the intensity ratio between peak and dip, and 
the spectral shape in the displayed energy range — are 
perfectly described by the model function. 

The most important feature of the spectrum in Fig. 0] 
is the drop of the intensity around 11 meV below unity, 
which is also lower than the BCS intensity and the nor- 
mal state DOS in this energy range. Choosing a simple 
Einstein- like model phonon at energy Hu>e, it was demon- 
strated [13 that the dip in the DOS appears approxi- 
mately at a binding energy slightly beyond Hlue + ^q ~ 
corresponding to the maximum position in 3mA. This is 
in accordance with our observation where the dip appears 
at an energy slightly higher than Hujd + ^o = 9 meV. In 
the BCS model, the superconducting gap width Ao at 
T = is related to the transition temperature T c by 
the dimensionless parameter 2Ao/ksT c — 3.50. This pa- 
rameter amounts to 4.497 for Pb (T c = 7.19 K), much 
larger than the BCS value; Al (Tn = 1.18 K) as the pro- 
totype BCS material has 3.535 This means, the gap 
width of Pb is larger than for a BCS like superconductor 
with the same T c . As a consequence, the characteristic 
weak photoemission structure from the thermally occu- 
pied singularity above Ep, which can be observed e.g. 
at T m 0.6T C for VaSi [J, is not resolvable because the 
Fermi-Dirac distribution suppresses the spectral intensity 
at energies farther away from the Fermi level. 

In conclusion, the high-resolution photoemission data 
presented in this paper clearly demonstrate the influence 
of the strong electron-phonon coupling on the electronic 
structure of lead in both the normal and superconduct- 
ing state. The experimental superconducting density 
of states, the renormalized band dispersion, the energy 
and the temperature dependence of the quasi-particle 
linewidth are in very good agreement with theoretical cal- 
culations based on the Eliashberg theory and represent 
consistently — and even quantitatively — the electron- 
phonon corrections in spectroscopic data. We hope that 
these results may also help to reveal the driving mecha- 
nisms behind the unconventional superconductors, like 
high-T c systems, superconducting heavy-fermion com- 
pounds and even superconducting organic materials. 
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